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Abstract 


The objective of this study is to investigate the transient behavior of a PEM fuel cell by using a one-dimensional, two-phase mathematical 
model. This model treats the catalyst layer as a spherical thin film-agglomerate. Effects of various transport parameter as well as other factors 
such as catalyst loading, gas diffusion layer thickness and liquid water permeability on the transient evolution of major model properties and 
cell performance are investigated in detail. Numerical results show that the evolution of ionic potential drop and oxygen consumption experience 
several steps before they reach steady state. The same situation can also be seen for the evolution of water saturation and current density. A detail 
inspection of these phenomena shows a close relation between these transport variables and cell performance. Parametric studies of other design 
factors’ effects reveal there optimum values which lead to a greater current output during its evolution period. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


It is of great important to increase the performance of fuel 
cells (FCs) for modern consumer application such as portable 
and mobile devices. A number of studies have shown that the 
cathode is the most important fuel cell component, and its flood- 
ing problem is the major cause of performance deterioration 
in proton exchange membrane (PEM) fuel cells. One way to 
alleviate this phenomenon is by developing a better water man- 
agement strategy. However, this requires a better understanding 
of electrochemical reactions and liquid water transport in the cat- 
alyst layers, which are the most complex components of the fuel 
cell. 

Due to the extreme small geometry of the electrodes, it is 
difficult to directly measure the liquid water distribution in the 
gas diffusion layers and catalyst layers. Many researchers have 
developed mathematical models to investigate how liquid water 
transport affects the PEM fuel cell performance. Springer et 
al. [1] first developed an isothermal, one-dimensional, steady- 
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state model for a complete polymer electrolyte fuel cell model. 
The model predicted an increase in membrane resistance with 
increased current density and demonstrated the great advantage 
of a thinner membrane in alleviating the resistance problem. 
Bernardi and Verbrugge [2,3] proposed a one-dimensional math- 
ematical model of the solid-polymer-electrolyte (SPE) fuel cell 
for liquid water transport in porous electrodes assuming a con- 
stant liquid water volume fraction and no interactions between 
liquid and gas flows. A numbers of catalyst layer models had 
been developed, include the interface model [1-3], the thin film 
model [4], the agglomerate model [5-11], and the thin film- 
agglomerate model [12-14]. Using the agglomerate model, the 
researchers at National Research Council of Canada studied 
the optimum performance of PEM fuel cell for a number of 
optimization parameters: type of agglomerate, catalyst layer 
thickness, catalyst layer porosity, distribution of Nafion con- 
tent, Pt loading, etc. [15-21]. Besides these models, the works 
of Wang and Wang [22,23] treated the catalyst layer as an indi- 
vidual zone with various conservation equations employed in 
the modeling of transient study and various time constants for 
the transient transport phenomena were proposed. 

As for the investigation of two-phase flow, Wang et al. [24] 
pioneered the research on this issue through analytic and numer- 
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concentration of species j (mol cm~?) 
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Subscripts and superscripts 
agg agglomerate 
CL catalyst layer 
effective 
gas phase 
gas diffusion layer 
membrane 
Nafion phase 
oxygen 
proton 
platinum 
saturation 
vapor water 
liquid water 


ical methods. A threshold current density was proposed to 
distinguish the scenarios between single- and two-phase regimes 
of water distribution and transport. In the subsequent works of 
their group, sophisticated models [25-28] were developed to 
simulate the flooding and liquid water distribution in PEM fuel 
cells. The mature multiphase mixture (M2) formulation with 


single set of conservation equations was employed to mimic 
the two-phase transport process. The results showed that it is 
important to consider the gas-phase pressure field in the GDL as 
a new mechanism of reactant transport. Another analytic model 
performed by Pasaogullari et al. [29], focusing on the liquid 
water transport in hydrophobic electrode, indicated that capillary 
transport dominates the water remove from flooded GDL. The 
model of Natarajan and Van Nguyen [30,31], which considered 
the evaporation or condensation of liquid water, also demon- 
strated the importance of its transport on cell performance. 
Moreover, comprehensive review summarizing the current stat- 
ues of the fuel cell modeling as well as point the development 
direction of this emerging field is presented by Wang [32]. 

With the consideration that the catalyst layer is the most 
important component in the electrochemical reaction of a fuel 
cell, also the experimental results of TEM images by Siegel et 
al. [7] clearly showed that carbon-supported, spherical Pt pellets 
exist in the catalyst layers. To mimic the actual morphology of 
the catalyst layers, a spherical thin film-agglomerate expression 
along with a one-dimensional, two-phase model is developed 
in this paper to investigate the transient evolutions of various 
model properties as well as cell performance. The boundary 
conditions at the CL/membrane interface also take into account 
the effects of water content on certain design parameters such as 
GDL permeability, CL permeability, catalyst loading, and GDL 
thickness. 


2. Mathematical modeling 


This study presents a transient analysis of a PEM fuel cell 
cathode. The model domain consists of a gas diffusion layer 
(GDL), a catalyst layer (CL), and a membrane, arranged as 
shown in Fig. 1. Our previously published cylindrical thin film- 
agglomerate model [13] is extended to a spherical model in 
this work. The electrochemical kinetics and transport of oxy- 
gen, vapor water, and liquid water are also analyzed in the GDL, 
CL, and proton exchange membrane (PEM). When the PEM 
fuel cell is turned on, air diffuses through the GDL from the 
channel to the catalyst layer. Meanwhile, oxygen dissolves into 
the Nafion film and eventually reaches the pellet surface, where 
liquid water is generated. The following assumptions are made: 


(1) The gas-phase obeys the ideal gas law. 

(2) Electronic resistance is negligible in the GDL and CL. 

(3) CL and GDL are assumed hydrophilic. 

(4) The catalyst pellets are spherical, consist of carbon- 
supported platinum and Nafion, and are covered by a Nafion 
film. 

(5) All the catalyst pellets in the CL have the same radius and 
Nafion film thickness. 

(6) The catalyst pellets are homogeneous, and oxygen diffuses 
into them via the Nafion film. 

(7) The oxygen reduction reaction generates only liquid-phase 
water. 

(8) The ionic potential between the anode catalyst layer and 
the membrane is approximately zero, because the hydrogen 
oxidation reaction is fast. 
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Fig. 1. Schematic of the model [12]. 


Our previous paper [13] provides details on the derivation 
of the model. The conservation equations applied are listed in 
Table 1, for the following five variables: 


(1) The concentration of oxygen in the gas-phase, [C6,]. 


(2) The concentration of water vapor in the gas-phase, [c$]. 
(3) The liquid water saturation level, [s] (the ratio between liq- 

uid water volume and void volume in the porous medium). 
(4) The concentration of liquid water in the Nafion phase, [C ]. 
(5) The ionic potential, [ġ] (in the Nafion phase). 


Experimental studies [5,7] have shown that the carbon- 
supported Pt catalysts in the CL are approximately spherical. 
As mentioned above, this investigation treats them as perfect 
spheres. The derivation of spherical thin film-agglomerate model 
described here is based on the cylindrical models of Lin et al. 
[12] and Chang and Chu [13]. The principal difference between 
the cylindrical model and the spherical model is the outer sur- 
face area per unit volume of the agglomerate. The surface area 
per unit volume of the agglomerate is defined as, 


Agglomerate’s surface area 


= 1 
i Geometric volume (1) 
For spherical thin film-agglomerate model, a, is: 
E 4n(rp + ôn)? 
(4/3)(rp + bn)? /(A — 5") 2) 
3 cL 
=- ae 
Cpt 


The initial values of all variables are zero, except for the liq- 
uid water concentration in the Nafion phase which is assumed 
to be in equilibrium with the water vapor activity in the gas- 
phase. All boundary conditions are listed in Table 2. The 
CL/membrane interface boundary condition for the concen- 
tration of liquid water in the Nafion phase (Eq. (3)) is taken 
from the model of Lin and Nguyen [14], where the water con- 
tent below 14 mol H20 (mol SO3~)~! was determined by water 
vapor activity. When the gas stream is saturated, however, and 
liquid water exists in the pores of the CL (between 14 and 
16.8 mol H20 (mol S037)7!), the water content varies linearly 
with the liquid water saturation level [1]. Since liquid water is 


generated when a fuel cell starts up, any liquid water in the 
CL/MEM interface pores can dissolve into the Nafion phase. 
This leads to the boundary condition: 


CNea — (0.043 + 17.81% — 39.84? + 36.007)C; + 2.8Cfs, 
(3) 


where a is the water activity in the gas-phase in the CL cath- 
ode and Cp is the concentration of fixed charged sites in the 
membrane. 
All the correlations used in this model are listed in Table 3. 
The governing equations and boundary conditions are dis- 
cretized by the finite difference method. The convergence 
criteria for a steady state are: 


yrew _ yold 
L 


za ei (4) 
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where V; is an arbitrary variable. 


3. Results and discussion 


Water management plays an important role on the perfor- 
mance of PEM fuel cells. Understanding their starting process 
and performance evolution are also crucial in the actual operat- 
ing applications. Simulation results for this issue by employing 
the PEM fuel cell cathode model is presented in the follow- 
ing sections. Investigations on how the structural parameters 
of the cell such as GDL permeability, CL permeability, cata- 
lyst loading, and GDL thickness affect its performance are also 
demonstrated. 

To validate the numerical model, the comparison of calcula- 
tion results and Navessin et al.’s experimental results are carried 
out. The comparison results are shown in Fig. 2. The experi- 
mental operating conditions and parameters of Navessin et al. 
are listed in Table 4, and those of the model are listed in Table 5. 
The /-V polarization curve represents the fuel cell’s steady state 
at an operating temperature of 25°C. The PEM fuel cell oper- 
ates under uniform conditions from the time it is turned on until 
a steady state is achieved. It is found that current model results 
agree well with the experimental data. 
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Fig. 2. A comparison of the model to experimental data. 


Fig. 3 displays the evolution of the total ionic potential drop 
(IPD) in the catalyst layer, for three different cell potentials Vs. 


x=CL 


IPD = Š. (plx=cr — $p). (5) 


x=GDL 


This plot can be divided into five sections: (1) a rapid rise, which 
lasts from startup to 0.01 s; (2) from 0.01 to 0.1 s, the potential 
plateaus (remains approximately constant); (3) at 0.1 s the poten- 
tial rises rapidly again to its maximum value; (4) the potential 
drops again; (5) a second plateau is reached representing the 
steady state. 

The total dimensionless oxygen consumption (TOC) in the 
catalyst layer also evolves over time, as shown in Fig. 4. This 
quantity is defined as: 


x=CL P8 g 
Co, |x=GDL — Co, 
TOC = ) call : (6) 
x=GDL O2 


The plot can be divided into six sections: (1) a rapid rise from 
startup to 5x 107? s; (2) the first plateau, which lasts from 
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Fig. 3. Evolution of the total ionic potential drop in the catalyst layer under 
various cell voltages. 
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Table 2 

Boundary conditions 

Variables X=0 GDL/CL CL/MEM X=L 

Cos c, = CBs No, lant = NOs lant No, lc. = 0 - 

c Cy = C" NŽlanL = N |epL Nle =0 - 

s s=0 Nwlapt = Nwlape Nwlet = 0 - 

ra = NA lao = 0 Cx“? = (0.043 + 17.81 — 39.807 + 36.003)Cr + 2.8Cfs a= as 
$ - iplapL = 0 iplcL = ipler g=0 

Table 3 


Correlations used in this study 


Diffusivity of oxygen in Nafion membrane [3] (cm? s7!) 
Diffusivity of oxygen in liquid water (cm? s7!) 
Diffusivity of oxygen in gas stream [33] (cm? s7!) 


Diffusivity of vapor water in gas stream [33] (cm? s7!) 


Henry’s constant for oxygen between Nafion and air [3] (cm? mol!) 


Henry’s constant for oxygen between liquid water and air [34] (cm? mol7!) 


Water content 
Net electro-osmotic drag coefficient [1] 
Conductivity of Nafion membrane [1] (1(Q cm)~!) 


Diffusivity of liquid water in Nafion membrane [1] (cm? s7!) 


5 x 107° to 5 x 107? s (consumption increases slightly during 
this phase); (3) another rise from 5 x 107? to 0.5 s; (4) a sec- 
ond plateau from 0.5 to 2s (where the consumption decreases 
slightly); (5) a rapid rise to the peak value from 2s onward; 
finally (6) a steady-state plateau. Oxygen consumption increases 
less than 15% in the first plateau, and decreases by only 2% in 
the second plateau. As the IPD increases, the fuel cell must over- 
come higher and higher activation energies to maintain a higher 
rate of electrochemical reaction. After the fuel cell overcomes 


Table 4 
Parameters used by Navessin et al. [15] 


Gas diffusion layer properties 
Porosity 
Thickness (cm) 


Catalyst layer properties 
Porosity 
Thickness (cm) 
Catalyst loading (mpt) (mg) 
Area per unit volume* (m7!) 
Volumetric fraction of Naion in catalyst pellet (R) 
Raduis of catalyst pellet, (Ragg) (cm) 
Exchange current density? (mA cm~?) 
Reference concentration of O2, (Co, ref)? (umol cm?) 


Membrane properties 
Porosity 
Thickness* (cm) 


Operation conditions 
Temperature (°C) 
Pressure (atm) 


* List value for HMEA-2 and HMEA-3, respectively. 


0.6 
0.035 


0.55 
0.005 
0.78 
727,664 
0.2 

2.5 x 1075 
1.27, 1.66 
2.34, 2.25 


0.35 
0.0062, 0.0061 


25.0 
1.0 


DN, = 0.0031 exp (— 2268) 
DW, = 2.14 x 107 “Gi. © 
DE = O1773( ae) 


Dy = 0. 256 ( sais, 15; a 
HÄ, = 1.33 x 10°exp ( =) 
HÄ, = 1.33 x 10%exp (=) 


cN 
La 
G 
na = 35, 
1 1 
kn = exp [1268 (355 — ne 0.005139. — 0.00326) 
Chasiz3 = inap — Ł)] 3.1 + 2.04) 


ana 
= es 


[241 
Craeaes) = 1076 exp [241 
cN 4 


Table 5 
Parameters used in this simulation (base case) 


Nasa) = 10~%exp [2416 ans — +) | (2.563 — 0.332 + 0.02464? — 0.000671) 


Gas diffusion layer properties 
Porosity 
Thickness (cm) 
Permeability (cm?) 
(-4 dpe. æ) (Dyne cm -2) 


Catalyst layer properties 
Porosity 
Thickness (cm) 
Permeability (cm?) 
(-4 dpe e) (Dyne cm -2) 
Ca loading (mp;) (mg) 


Specific surface area of Pt (apt) (cm?(mg P!) 
Volumetric fraction of Nafion in catalyst pellet (G9) 


Radius of catalyst pellet (Ragg) (cm) 
Thickness of Nafion (cm) 
Exchange current density (io,ref 0°C) (A cm~?) 


Membrane properties 
Thickness (cm) 
Porosity 
Fixed charge site concentration (mol cm~’) 


Operation conditions 
Temperature (°C) 
Pressure (atm) 
Mole fraction of O2 in inlet 
Humidity in inlet (%) 
Humidity at anode (%) 
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Fig. 4. Evolution of the total dimensionless oxygen consumption in the catalyst 
layer under various cell voltages. 


the maximum activation energy, the IPD decreases. As shown 
in Fig. 4, when V; = 0.2 V and 0.4 V, after 2 s, the TOC increases 
rapidly. Thus, after a fuel cell overcomes the maximum activa- 
tion energy, it retains a high electrochemical reaction rate with 
a lower IPD. 

Fig. 5 shows the evolution of the liquid water saturation level, 
measured at various cell voltages. The liquid water saturation 
level plot can be divided into four sections. From startup to 
107°? s, liquid water accumulates in the catalyst layer. Between 
107? and 0.5s, liquid water begins diffusing into the gas dif- 
fusion layer. During the next phase, from 0.5 to 2s, the liquid 
water saturation level in the CL decreases as more and more 
liquid water occupies pores in the GDL. In the final phase, the 
liquid water saturation level increase in both CL and GDL until 
a steady state is reached. Although the liquid water saturation 
level decreases in the CL between 0.5 and 2s, note that it will 
increase again once the electrochemical reaction rate generates 
enough liquid water to fill the Nafion phase. When V, =0.6 V, 
for example, the electrochemical reaction rate is moderate and 
less liquid water is generated. In this case, at the time after 2 s, 
the liquid water saturation level decreases (Fig. 5c). In contrast, 
the liquid water saturation level increases after 2s for V; =0.2 
and 0.4 V (Fig. 5a and b). 

The fuel cell cannot be turned on if the Nafion phase is com- 
pletely dry, so it is assumed that the concentration of liquid water 
in the Nafion phase is equal to the concentration of fixed charge 
sites at the anode. Fig. 6 shows the evolution of the CL water 
content over time at three different cell voltages. These plots can 
be divided into three sections: from startup to 3 s, the water con- 
tent decreases; from 3 to 13 s, the water content increases; and 
in the final phase, the water content is constant. Initially there 
is less liquid water amount in the gas pores than in the Nafion 
phase, so at first liquid water flows from the Nafion phase to 
the gas pores. When more liquid water begins to be generated, 
it dissolves into the Nafion phase again and we see increasing 
water content in both the CL and the membrane. 

The current density evolution for various cell voltages is 
shown in Fig. 7. These plots can be divided into five sections. 
From startup to 5 x 1073s, the current density rises rapidly. 


T 


ooooo0oo00o0 
NOoOhWND— OO 


0.44 


Liquid water saturation level, s 


0.1 


0.0 


o 
a 


gooocosso 
BnRwonr Ao 


o 
P 


Liquid water saturation level, s 
oO O 
N w 


9 
A 


0.20 } 


oO 
=à 
ao 


Liquid water saturation level, s 
o 
Td 
Oo 


o 
o 
a 


0.00 


30 
25 20 15 49 n 
Time [sec] P 


Fig. 5. Evolution of the liquid water saturation level under various cell voltages, 
(a) Vs =0.2 V, (b) V; =0.4 V and (c) V; =0.6 V. 
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Fig. 6. Evolution of the water content under various cell voltages, (a) Vs = 0.2 V, 
(b) V; =0.4 V and (c) V; =0.6 V. 
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Fig. 7. Evolution of the current density under various cell voltages. 


From 5 x 1073 to 2 x 107° s, the current density plateaus. The 
current density begins to rise rapidly again at 2 x 107? s, even- 
tually reaching its maximum value, then drops briefly. Finally, it 
reaches a steady-state plateau. As seen in Fig. 5, liquid water hin- 
ders oxygen transport as it diffuses from the CL to the GDL. The 
electrical conductivity is also reduced during this period, due to 
decreased water content, and the current density levels off as a 
result. Note that the total oxygen consumption increases rapidly 
after 0.05 s (Fig. 4), and that the ionic potential drop increases 
rapidly after 0.1s (Fig. 3). Thus, the current density increases 
rapidly after 0.1s. Oxygen transport resistance increases after 
0.1 s, however, because at this point the liquid water saturation 
level is increasing in both GDL and CL (Fig. 5). Thus, the current 
density levels off after 1 s. 

The current density drops off after reaching its maximum 
value for the following reason. The higher the current density, the 
more liquid water is generated and the more gas pores are occu- 
pied by liquid water. When this happens, however, the catalyst 
pellets also have less oxygen on their exposed surface. This tends 
to reduce the electrochemical reaction rate and causes a drop in 
the current density. At lower cell voltages the electrochemical 
reaction rate is higher, so more liquid water is generated and 
we see a more extreme drop in the current density. For instance, 
when V, =0.2 V the drop is more than 30%. 

Fig. 8 shows the effect of modifying GDL permeability for 
V, =0.4 V. Cases 1 through 4 used permeabilities of 9 x 10719, 
1 x 107°, 1 x 1078, and 1 x 1077 cm?, respectively. The higher 
the permeability, the more liquid water will be drained; this low- 
ers the amount of liquid waver available and reduces the current 
density. If the GDL permeability is too low, however, liquid 
water cannot drain effectively. This will reduce the available 
oxygen on the catalyst surface, decreasing cell performance and 
limiting mass transport. 

Fig. 9 examines the effect of CL permeability when 
V; =0.4 V; cases 1 through 4 use values of 8 x 107!2,3 x 107!!, 
3 x 107!°, and 3 x 107° cm?, respectively. Cases 3 and 4 have 
the highest permeabilities and effective water drainage. But 
when liquid water drains from the fuel cell quickly, it also 
reduces the water content of the Nafion phase. Thus, a lower 


S.-M. Chang, H.-S. Chu / Journal of Power Sources 172 (2007) 790-798 797 


Current density [A cm-2] 


“104 103 10-2 10 100 101 102 108 
Time [sec] 


Fig. 8. The effect of GDL permeability on the current density (V; =0.4 V). 


CL permeability means that there will be more liquid water in 
the Nafion phase. The maximum current densities achieved in 
cases | and 2 are similar, but the steady-state current density 
in case 2 is higher. This is because in case 1, liquid water is 
occupying more of the gas pores and effectively covering the 
catalyst pellet surface, causing a series concentration overpo- 
tential. 

The current density evolution for V;=0.4V is shown in 
Fig. 10, under catalyst loadings of 0.2, 0.3, 0.4, and 0.5 mg for 
cases | through 4, respectively. A higher catalyst loading means 
a more rapid electrochemical reaction. The current densities in 
cases 1 and 2 are obviously lower than those in cases 3 and 4, 
because the catalyst loading is too low in the former to drive the 
reaction. The maximum current density in case 4 drops rapidly, 
on the other hand, because of fuel cell flooding. Oxygen cannot 
reach the surface of the catalyst pellets when there is too much 
liquid water. This suggests that excessive catalyst loading does 
not improve cell performance. 

Finally, Fig. 11 demonstrates the evolution of current den- 
sity for Vs =0.4 V under various thicknesses of the gas diffusion 
layer. Through case | to 6, the thicknesses are 220, 235, 240, 245, 
250, and 265 um, respectively. The current density increases to a 
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Fig. 9. The effect of CL permeability on the current density (V; =0.4 V). 
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Fig. 10. The effect of catalyst loading on the current density (V; = 0.4 V). 
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Fig. 11. The effect of GDL thickness on the current density (V, = 0.4 V). 


constant (steady state) value after 2 s for GDL layers less than or 
equal to 240 pm thick. For the remaining cases, the current den- 
sity rises again after 2 s, reaches its maximum, then decreases 
in value and reaches a steady state. A thicker gas diffusion layer 
means a higher maximum current density, but this value may 
overshoot the sustainable current. Thicker GDL also imposes 
a greater resistance to liquid water transport, which promotes 
flooding of the fuel cell. Liquid water drains more quickly with 
a thinner gas diffusion layer, but the fuel cell cannot be too dry 
or its performance will deteriorate. 


4. Conclusion 


A transient, one-dimensional, two-phase model of a PEM 
fuel cell cathode is presented in this study with the catalyst layer 
treated as a spherical thin film-agglomerate. This model is used 
to investigate the transient transport of gaseous species, protons, 
and liquid water. The effects of GDL permeability, CL perme- 
ability, catalyst loading, and GDL thickness on cell performance 
are also investigated. 

Our model results agree well with the experimental data. The 
observed IPD and TOC variations show that after the fuel cell 
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overcomes its maximum activation energy, a high electrochem- 
ical reaction rate with low IPD is maintained. Electrochemical 
reactions affect the amount of liquid water in the cell; when less 
water is generated, the cell reaction rate is reduced. Neverthe- 
less, the higher liquid water content and saturation level tend to 
lead to the flooding catastrophe. 

Parametric sensitivity analysis also shows that at higher per- 
meability, more liquid water is able to drain from the cell. Cell 
performance will decrease if there is too little liquid water, but 
may also be influenced if the water cannot drain effectively. 
Among the values adopted in this paper, the optimum GDL 
permeability is Keo. = 1 x 107° cm? and the optimum CL 
permeability is KS = 3 x 1071! cm’. 

Moreover, a cell with higher catalyst loading leads to more 
rapid electrochemical reactions, but this may also result in too 
much liquid water generation. This ultimately impedes oxygen 
transport, so excessive catalyst loading does not improve cell 
performance. Among the values considered in this study, the 
optimum loading is mp; = 0.4 mg. 

Finally, GDL thickness was found to play an important role 
in liquid water transport. A thinner layer offers less resistance 
to liquid water diffusion and drainage. Water content and cell 
performance therefore tend to remain constant. When the gas 
diffusion layer is thick, liquid water cannot drain effectively and 
the catalyst pellets will have less surface oxygen. According to 
these results, the optimum thickness is dgpr = 245 um. 

Fuel cell performance and durability are also strongly influ- 
enced by impurities in the hydrogen gas used. A particular level 
of carbon monoxide yields a stable performance loss. There- 
fore, another research topic which is currently undergoing in 
our group is the investigation of the CO poisoning effect with 
the thin film-agglomerate model. 
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